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In this paper, we compare dry-snow extinction coefficients 
derived from radar altimeter data with brightness temperature 
data from passive microwave measurements over a portion of 
the East Antarctic plateau. The comparison between the 
extinction coefficients and the brightness temperatures shows 
a strong negative correlation, where the correlation coefficients 
ranged from -0.87 to -0.95. The extinction coefficient of the 
dry polar snow decreases with increasing surface elevation, 
while the average brightness temperature increases with surface 
elevation. Our analysis shows that the observed trends are 
related to geographic variations in scattering coefficient of 
snow, which in turn are controlled by variations in surface 
temperature and snow accumulation rate. By combining 
information present in the extinction coefficient and brightness 
temperature datasets, we develop a model that can be used to 
obtain quantitative estimates of the accumulation rate of dry 
polar snow. 

INTRODUCTION 

The role of the ice sheets in global processes has become a 
major focus of research within the last two decades. Climatic 
warming caused by an increase in levels of carbon dioxide and 
other greenhouse gases has become a major political and 
environmental issue. Although the ice sheets play a critical 
role in global climatic processes, they remain one of the most 
poorly explored regions on the Earth. Satellite remote sensing 
systems provide the only proven means for observing the polar 
regions with the temporal resolution and spatial coverage 
required for meaningful ice-sheet studies [1]. Thus, it is 
important to compile new time series from existing satellite 
data and to continue the existing time series of data. This will 
provide a baseline set of geophysical measurements to improve 
our understanding of the role of polar processes in global 
climate. The baseline time series can then be analyzed to 
determine variations in geophysical parameters as evidence of 
climatic anomalies, trends, and even climatic change. 

Passive microwave radiometers have been used to produce 
images of the polar regions since the early 1970s. Spatially 
varying patterns of microwave brightness temperature have 
resembled geographic variations in snow accumulation in some 
areas and physical temperature in others. For dry polar snow, 
the microwave emissivity depends strongly upon the snow grain 
size distribution with depth, which is determined primarily by 
accumulation rates and snow temperature. Using radiative 
transfer theory, Zwally [2] and Comiso et al. [3] correlated 
snow grain size with measured emissivity to calculate accumu- 
lation rates at a few selected sites in Antarctica. Tiese 
investigations showed that emitted energy from depths up to 20 
m contribute to the observed radiometric brightness tempera- 
ture. 

Satellite radar altimeters have been used to measure surface 
elevations over the ice sheets. Recent work has shown that tlie 


shape of altimeter return waveforms from the ice sheet is 
significantly affected by sub-surface volume scattering [4], [5]. 
By using a combined surface and volume-scattering model, 
Davis and Zwally [6] extracted dry-snow extinction coefficients 
from radar altimeter data. The extinction coefficient of snow 
characterizes the amount of signal penetration beneath the ice- 
sheet surface and is determined by its absorption and scattering 
properties. The smaller the extinction coefficient, the greater 
the amount of altimeter signal penetration. 

SATELLITE DATA COMPARISON 

In this section we compare active and passive microwave 
measurements over the Antarctic continent and show that 
geophysical properties of the ice sheet can be studied using 
these measurements. For this study, we selected a region from 
the East Antarctic ice sheet that covered a wide range of 
surface elevations. Datasets from the Geosat radar altimeter and 
the Scanning Multichannel Microwave Radiometer (SMMR) 
were obtained for the summer and winter months of 1986 and 
1987 (Table I). The datasets are in a region bounded by the 
latitudes 71.5 to 72.1°S and the longitudes 80 to 100°E. An 
elevation profile for the selected region is shown in Figure 1. 
The profile begins at an elevation of 2550 m at 80°E, increases 
to an elevation of 3150 m at the ice divide (94°E), and then 
decreases to an elevation of 3000 m at 100°E. 

Following the approach by Davis and Zwally [6], we fit a 
combined surface and volume-scattering model [5] to Geosat 
altimeter return waveforms in 1° subdivided regions for the 
two-month periods corresponding to the summer and winter 
seasons. After fitting the individual waveforms, we averaged 
the extinction coefficient (k^ values from the model fits 
together for each subdivided region. By fitting large numbers 
of individual waveforms, sufficient samples of the extinction 
coefficient are obtained to characterize seasonal and geographic 
differences in the surface properties of the ice sheet. For the 
SMMR satellite, we selected 10-GHz data since they are closest 
to the 13.5 GHz operating frequency of the Geosat altimeter. 
In this region the SMMR satellite made brightness temperature 
measurements every other day. Measurements at a selected 
longitude in each 1° region were averaged together over two- 
month periods to obtain the corresponding seasonal average. 

A plot of the extinction coefficient vs. longitude for the 
summer and winter seasons of 1986 is shown in Figure 2. The 
data point at each longitude represents the sample mean for the 
region extending from the given longitude to a longitude one 
degree higher; i.e., the data point at 89°E corresponds to the 
region from 89-90°E. The extinction coefficient decreases from 
a maximum value of 0.2 m' 1 at 80°E to a minimum value of 
approximately 0.08 m" 1 at 99°E This corresponds to penetra- 
tion depths that range from 5.0 to 12.5 m. The 1987 data show 
the same extinction coefficient variation with longitude. The 
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Table l. Geosat and SMMR Datasets 

Area Years 

Summer 

Months 

Winter 

Months Data Region 

East Antarctica 1986,1987 

Jan -Feb 

July- Aug 71.5-72.1* S. 80-100° E 9 


* arc-like region bounded in latitude 



Figure 1. Elevation profile for the data region in East Antarcti- 
ca from 80 to 100°E at a latitude of 7 2°S. 


decrease in k e from 80-93°E is strongly related to the increase 
in elevation from 2550-3125 m (Figure 1). Beyond the ice 
divide at 94°E, the extinction coefficient experiences a sharp 
decrease even though the elevation drops slightly from 3150- 
3000 m. This decrease is probably related to a marked 
difference in the physical properties of the snow on either side 
of the ice divide. 

The large-scale trend in the extinction coefficient is strongly 
related to surface elevation. Surface elevation can be used as 
a proxy for studying the effect of the mean annual surface 
temperature on the geophysical properties of the snow. As 
elevation increases, the mean annual surface temperature 
decreases. The growth rate of ice crystals in polar snow 
depends upon the mean annual surface temperature, where 
warmer temperatures produce higher growth rates. Larger 
growth rates increase the mean size of the snow grains, which 
increase the extinction coefficient. Thus, regions that have 
lower surface elevations generally have higher extinction 
coefficients. 



the same increase with longitude. Theoretical modeling has 
shown that the microwave emissivity is controlled by variations 
in the density and grain size of polar fim as a function of depth 
[2], [3]. As the average grain size increases, the amount of 
emitted energy from greater snow depths is reduced because the 
upwelling radiation is attenuated due to an increase in the 
scattering loss. This causes the observed brightness tempera- 
ture to decrease. Thus, the increasing brightness temperature 
with longitude suggests that the average grain size is decreas- 
ing. 



Figure 3. Mean brightness temperature vs. longitude for the 
summer (boxes) and winter (triangles) seasons of 1986. 

Comparing the brightness temperature and extinction coefficient 
plots, we see that as the extinction coefficient decreases, the 
brightness temperature increases. Smaller average grain sizes 
decrease the extinction (primarily scattering) coefficient of 
snow. Thus, emitted energy from greater snow depths increas- 
es, thereby increasing the observed brightness temperatures. 
Linear correlation coefficients, r, were computed for the 
extinction coefficient and brightness temperature variations and 
these are summarized in Table II. The correlation coefficient 
varies from -0.87 to -0.95, indicating a strong negative correla- 
tion. The negative correlation is due to the opposite effect that 
the snow grain size has on the observed brightness temperatures 
and extinction coefficients. 


Table It. - k e Correlation Coefficients 

Summer *86 Winter *86 

Summer ’87 

Winter *87 

Correlation Coefficient -0.92 -0.94 

-0.87 

-0.95 


While the large-scale trends in the extinction coefficient and 
brightness temperature datasets are clearly dominated by the 
snow grain-size dependence upon mean annual temperature, and 
hence elevation, there are some regional deviations from the 
overall trend. Note that in the region from 88-93°E, both the 
extinction coefficient and mean brightness temperature are 
nearly constant as the elevation increases from 2950 to 3100 m 
for both the summer and winter seasons. Although the 
correlation coefficients computed for the extinction coefficient 
and brightness temperature variations with elevation are large, 
there are significantly smaller ( I r I « 0.76 to 0.84) than those 
in Table II. This suggests that the extinction coefficient and 
brightness temperature variations are influenced by another 
geophysical parameter. 


Figure 2. Extinction coefficient vs. longitude for the summer 
(boxes) and winter (triangles) seasons of 1986. 

A plot of the 1986 mean brightness temperature vs. longitude 
is shown in Figure 3. For both seasons, the mean brightness 
temperature increases with longitude. The 1987 data showed 


SYNTHESIS OF GEOSAT ALTIMETER AND SMMR DATA 

Work by Gow [7], [8] showed that the growth of snow grains 
in polar snow is not only dependent upon mean annual temper- 
ature, but also on the accumulation rate. For small accumula- 
tion rates, the snow grains are exposed to larger surface 
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temperatures for a longer period of time. For high accumula- 
tion rates, the snow grains are buried at a greater rate and are 
insulated from the large temperature gradients at the surface. 
The net result is that the recrystallization rates of the snow 
grains increase as the accumulation rate decreases. Based upon 
this work, Zwally [2] showed that the equation 

r (z) 3 = r 0 3 + k o ^e x p[^§] (1) 

accurately fit crystal size data from pit studies, where r(z) is the 
crystal radius at depth z, Tq is the radius at the surface, q is the 
mean snow density, A is the mean annual accumulation rate, D 
is the activation energy, R is the gas constant, and T is the 
mean annual surface temperature. Here we see that as T 
increases, the growth rate increases, and conversely, as A 
increases, the growth rate decreases, just as described above. 
The extinction coefficient of snow is a linear combination of 
the absorption and scattering coefficients 

ke = k a + k s , (2) 


where the scattering coefficient, k^, in turn depends upon the 
radius of the snow grain. A snow grain is usually comprised 
of one or two ice crystals and Gow [7] showed that the ratio is 
approximately 1.4:1. The absorption coefficient depends upon 
frequency, snow density, and the imaginary part of the dielec- 
tric constant of ice. Using dielectric properties reportal by 
Matzler [9] and assuming a mean density of 0.45 g/ cm , we 
calculated an absorption coefficient value of k t = 0.04 m . 
The mean snow density over this region is nearly constant, and 
so the absorption coefficient is also assumed to be constant. 


The brightness temperature measured at the surface of the ice 
sheet, T b , is found by solving the radiative transfer equation 
and is given by 


T b = (l-R) secG J [ T(z) k a (z) +T sc (z) k s (z) ] 1 


(3) 


exp |-^9/(k a (z>k s (z’))dzj dz 


where R is the power reflection coefficient, 0 is the angle of 
transmission, T(z) is the physical temperature at depth z, and 
T (z) is the scattering source term at depth z. This equation is 
extremely complicated to solve and requires detailed knowledge 
of the depth variations of temperature, the scattering source 
term, and the absorption and scattering coefficients. Previous 
investigations have used complicated numerical methods to 
solve this equation and relate the observed brightness tempera- 
tures at the ice-sheet surface to the underlying geophysical 
properties. 

Here we propose a simple technique which can be used to 
relate T b to snow grain size without having to numerically 
solve (3). The technique is based upon the fact that the 
altimeter data provide an independent measurement of the 
extinction coefficient k e . The depth of penetration, dp, i$ 
inversely related to the extinction coefficient by 


This provides a reasonable estimate of the depth over which the 
various parameters in (3) contribute to the observed brightness 
temperature. We can solve (3) by assuming that all the terms 
inside the integral are constant with depth to obtain 


T b = (l-R) 


' *<3 - 

, k a +k s T + 




(5) 


Very little is known about the scattering source term, T^, but 
using the extinction coefficient data from the altimeter measure- 
ments we can approximate T sc using 

Tsc = TV^[t-T 5 ] . (6) 

As k approaches k a (no scattering), T^ becomes 0, and as 1^ 
approaches 0, T b approaches T^, as it should. In developing 
(5), it is assumed that the temperature is constant with depth, 
which is an approximation to the physical situation. Therefore, 
instead of using a surface temperature in (6), we calculate an 
effective mean snow temperature <T> by integrating a tempera- 
ture profile model used by Zwally [2] over the depth of 
penetration given in (4), i.e. 

d P 

<T> = / T(z) dz . (7) 

0 

This is then used along with the extinction coefficient data to 
estimate the scattering source term in (6). Then, we can solve 
(5) for k,. relating it to the brightness temperature, T b . Using 
the same form as (7), (1) is integrated to obtain the effective 
mean grain size, <r>, which is related to the mean surface 
temperature and accumulation rate. Assuming Rayleigh 
scattering and multiplying it by a dense medium coefficient f d 
= 0.3 [3], we can solve the integrated form of (1) for the mean 
accumulation rate. 

Figure 4 shows a plot of the mean accumulation rate vs. 
longitude. The accumulation rates vary from 4 to 18 g/ cm and 
agree with published maps of accumulation rates for the East 
Antarctic plateau. These maps are based upon sparse accumu- 
lation measurements from station data and show very little 
detail on a regional scale. The derived rates in Figure 4 show 
an excellent regional detail. From 80 to 93°E, we see a steady 
decrease in accumulation rate from 9 to 4 g/cm 3 . Also note 
that the minimum rate occurs at the ice divide (Figure 1). In 
this region, the opposing effects of mean annual temperature 
and accumulation rate compete for control of the mean grain 
size. The decrease in surface temperature with elevation causes 
the average grain size to initially decrease, which in turn causes 
the extinction coefficient to decrease and the brightness 
temperature to increase from 80 to 88°E, However, the 
decreasing accumulation rate with increasing elevation in this 
region competes with the temperature effect, where the lower 
accumulations tend to produce a larger mean grain size. The 
accumulation effect on grain size balances out the temperature 
effect on grain size in the region from 88 to 93°E, and we see 
a leveling off of the extinction coefficient and brightness 
temperature data. In the region east of the ice divide (> 93°E), 
we see a dramatic increase in the accumulation rate where the 
maximum value is four times the minimum value at the ice 
divide. In this region, the accumulation rate dominates the 
mean grain size. If temperature were the controlling factor, we 
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[ “ would expect to see a increase in the extinction coefficient and 

j a decrease in T b as the elevation decreases from 3150 to 3000 

» m. However, the exact opposite is true, indicating that the 

• higher accumulation rates on the east side of the ice divide 

j cause the mean grain size to continue to decrease, despite the 

* decrease in elevation. This accounts for the further increase in 

: brightness temperature and decrease in the extinction coefficient 

r (Figures 2 and 3). 
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Longitude 

ca 

Figure 4. Accumulation rate vs. longitude derived from 
combining the extinction coefficient and brightness temperature 
data. 

CONCLUSION 

Clearly, for the region of East Antarctica studied here, the 
microwave brightness temperature is strongly related to the dry- 
snow extinction coefficient. The synthesis of extinction 
coefficient and brightness temperatures data enabled reasonable 
estimates of the snow accumulation rate to be obtained. These 
estimates are the first to provide a detailed picture of small 
regional variations in the accumulation rate. The continuous 
time series of passive-microwave data (1978 to present) 
provides the opportunity to study long term variations in the 
geophysical properties of the continental ice sheets. Our ability 
to interpret these measurements can be significantly enhanced 
by including measurements from other satellite remote sensing 


systems. This type of data synthesis should be used to 

establish a baseline set of geophysical parameters estimates. 

This baseline set of measurements would be a valuable tool for 

studying long-term climatic variations of the ice sheets. 
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